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ABSTRACT 


The study of a biokinetic model development for bioremediation of 1kg of crude oil contaminated soil inoculated with crushed sun- 
dried and room-dried Azadirachta indica stem was conducted. First, physico-chemical analysis was carried out on the samples used 
in the experiment. The experiment was carried out ex-situ in twenty-six batch reactors, two which were the control. The other 
twenty-four reactors containing 1kg each of two types of soil: sandy in twelve reactors and loamy in another twelve were polluted 
with 100ml of Bonny light crude oil and inoculated with two samples of crushed Azadirachta indica stem. The quantities of 
Azadirachta indica stem were varied from 50g to 100g at 10g difference. The set-up was observed for 35 days as GC-MS analysis 
was done on the soil samples on weekly basis to determine crude oil depletion in each reactor. It was observed that contaminated 
soils treated with 100g blended room-dried Azadirachta indica stems had the highest TPH decline from 48508.92ppm to 
3707.78ppm for loamy soils (LRD) and 35818.69ppm to 1003.18ppm for sandy soils (SRD). Based on this observation, a predictive 
model was developed for either soil samples using both Analytical methods and numerical methods that predicted with high 
precision, the TPH content of the soil at any given time.TPH at 21 days was calculated to be 9085.47mg/kg which is similar to that 
obtained from experiment observed to be 9085.48mg/kg. The highest error was observed in day 7 where calculated TPH was 
23477.06mg/kg while the experimental value was 26865.7mg/kg. For LRD, TPH at 35 days was calculated to be 3707.81mg/kg which 


disc@@very 
‘Cotee the encudered 


Based 50 


ARTICLE 


is similar to that obtained from experiment observed to be 3707.78mg/kg. The highest error was observed in day 7 where calculated 
TPH was 35707.07mg/kg while the experimental value was 37851.62mg/kg. 


Key words: Biokinetics, crude oil, remediation, Dogoyaro, Azadirachta Indica, stem. 


1. INTRODUCTION 


There is no gain saying that the upstream, midstream and downstream operations of the oil sector is the paramount cause of 
environmental pollution in Nigeria. After oil discovery in 15" January 1956, there have been numerous oil-spill cases: between 1976 
and 1996 alone, a total of 4647 incidents resulted in the spilling of about 2,369,470 barrels of oil into the environment, of which an 
estimated 1,820,410.5 barrels (77%) were not recovered [1]. This amount of crude oil left unmitigated in the soil could pose very 
precarious threats to the environment and must be taken care of immediately in order to preserve aquatic or plant life and soil 
nutrients [2-3]. 

Paradigms for treating crude oil contaminated soils include incineration, soil vapor extraction, solidification/stabilization, soil 
washing, bioremediation etc. Some of these methods have their disadvantages because in the aftermath of the treatment, toxins 
noxious to plant and animal life are released into the soil [6-8]. Aside this, these methods are relatively expensive hence the need for 
cheaper alternatives, cost effective and very effective means of control [4-5]. In recent times, bioremediation has become one of the 
most propitious technologies with increasing demand for remediating crude oil contaminated soils because contaminants can be 
removed using innocuous and soil-friendly microbial colonies [6]. 

Numerous studies have established bioremediation as a cheap, eco-friendly and effective way of amending soils contaminated 
with crude oil. Studied conducted on plant extract of various species and the use of two species of bitter-leaf; Vernonia Galamensis 
and Vernonia Amygdalina, in the bioremediation of crude oil contaminated soils reveals their significant on crude oil remediation [7- 
8]. The leaves of each species were applied to the contaminated soil either room dried, sun dried or wet and blended [9]. The 
microbes (Pseudomonas Aeruginosa, Staphylococcus Aureus and Eureschia Coli) and phytochemicals in the leaf extracts were the 
agents responsible for the degradation of the metals as well as the hydrocarbon contents in the soil [10-12]. Results showed that the 
use of forty-gram of both species of Vernonia extracts reduced more than fifty-percent of the contaminants concentration in the 
soil, after forty days of investigation [13-15]. Therefore, the extracts of both species of Vernonia were effective as agents of 
bioremediation in a crude oil contaminated soil. However, studies conducted by various research groups recommended that, for the 
remediating effect of the either of the extract of both Vernoniaspecies (depending on which is available) to be maximized, it should 
be applied to the polluted soil wet and blended [16-20]. This is because the microorganisms present in the leaves to carry out the 
bioremediation are still active [21-24]. The result from this study showed that wet blended Vernonia extracts remediated more than 
fifty percent of the initial crude oil content of the contaminated soil [25]. 

Investigation on the possibility of enhancing crude oil bioremediation by supplementing soil with cost-effective organic 
materials derived from two widespread trees grown in Kuwait, Conocapus and Tamarix [25-26].The leaves of both plants were 
sterilized and blended and phytochemical analysis was carried out on the soil. Amendments in soils increased the soil’s microbiota 
count by up to Ninety-Eight percent and enhanced their activity by up to Ninety-Five percent [27-30]. 

Studieson the effectiveness of grounded neem tree leaves in crude oil bioremediation process. Results showed that the leaves of 
Azadirachta indica was able to support and promote the growth of hydrocarbon degraders (Pseudomonas species) indicating that 
they play complementary role in bioremediation [31]. The biostimulation potential of Azadirachta indica leaves in increasing the 
microbial population in crude oil contaminated vessel may be due to the nitrogen and phosphorus content of the stem [32]. 
Increase remegiation of polluted site was observed when using NPK fertilizer, urea fertilizer and poultry droppings effectively 
stimulated the growth of bacteria into utilization of crude oil [33].The research carried out established that the leaves of Azadirachta 
indica is very effective in bioremediation of crude oil contaminated soil by stimulating the growth of crude oil utilizing endophytes 
[34]. 

Obviously, in recent time enormous research had been carried out on bioremediation. However, not much work has been done 
in the mathematical modelling of the process involved in it. This study adopted a Monod-based-mathematical model to describe 
the bioremediation process of crude oil contaminated soils with Azadirachta indica stem as remedant [35-36]. 

Owing to the discrepancies in the physico-chemical properties of different soils and the difference in chemical constituents of 
different crude oil samples in different environments, previous models may not be adequate enough to describe the treatment of oil 
spill in the Niger Delta 37-38]. For Instance, models developed by some research groups, describes crude oil bioremediation in arid 
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Tehran may not be useful in describing crude oil remediation in Niger-Delta [39].This study aimed at designing a simple andrealistic 
model that describes the treatment of crude oil contaminated soils. 


2. MATERIALS AND METHODS 


Materials and Equipment Used for the Experiment 
The soil samples and the Azadirachta indica leaves were obtained from a farm within Rivers State University vicinity. 


Sample Preparation and Experimental Setup 

Twenty-six (26) reactors were used in the experiment of which 2 were control experiment containing crude oil contaminated soils 
without bioremediant. The other 24 reactors contained crude oil contaminated soils thoroughly mixed with varying quantities of 
either sun-dried or room-dried blended Azadirachta indica stem. The crude oil sample obtained was investigated for Total 
Petroleum Hydrocarbon content. The sandy and loamy soils were weighed 1kg into thirteen places for each type of soil. Some of the 
Azadirachta indica stem were sun-dried for 7 days and the others were room-dried for 14 days. Both the sun-dried and room-dried 
leaves were thoroughly blended.The weighed soil samples were contained in 26 reactors, 1kg of soil to 1 reactor. Each 1kg of soil 
sample was then contaminated with 100ml of crude by thorough mixing it with the soil. The blended sun-dried stems were weighed 
ranging from 50g to 100g with one replicate for each measure. The room-dried stems were weighed ranging from 50g to 100g with 
one replicate for each measure also. This summed up to 24 weighed blended stem samples. The weighed stemsamples were 
systematically applied on 24 of the contaminated soil samples. Below is a schematic representation of the experimental set up: 


Table 1: Schematic Representation of the Experimental Set-Up. 


50g 60g 70g 80g 90g 100g 
Sandy Loamy sandy loamy Sandy Loamy sandy loamy sandy Loamy sandy | loamy 
SSD LSD SSD LSD SSD LSD SSD LSD SSD LSD SSD LSD 
SRD LRD SRD LRD SRD LRD SRD LRD SRD LRD SRD LRD 


SSD represents the reactor containing sandy soil contaminated with crude oil and inoculated with blended sun-dried Azadirachta 
indica stem. SRD is the same with SSD except that blended room-dried Azadirachta indica stems were used instead of the sun-dried. 
To each type of soil contaminated with crude oil, quantities of either sun-dried or room-dried blended Azadirachta indica stem 
ranging from 50g to 100g were added as illustrated in table 1. Each combination was thoroughly mixed to ensure uniformity of the 
components. 

Before inoculation, the contaminated soils were analyzed for initial conditions like Total Hydrocarbon Content, initial microbial 
activity, conductivity, pH etc. Subsequently, the 24 experimental samples were analyzed every week for a period of 8 weeks. Within 
this period, 10ml of water was added to each of the 24 inoculated soil samples to aid easy movement of microbes. 


Gas Chromatography Procedure 

The soil sample was poured into a 1 liter separation funnel. 50mL of methylene chloride was added to the sample bottle seal and 
shook for 30 seconds to rinse the inner surface. The solvent was transferred to the separation funnel and then sample was extracted 
by shaking the funnel for 2 minutes with periodic venting to release excess pressure. 

The organic layer was allowed to separate from the water phase for a minimum of 10 minutes and then the methylene chloride 
was extracted in a 250mL flask. A second 60mL volume of methylene chloride was added to the sample bottle. The separation funnel 
and the column was rinsed with 20mL of the solvent into the extract. The extraction procedure was repeated a second time and 
combined with the other extract in an Erlenmeyer flask. The third extraction was performed in the same way. The combined extracts 
were poured through a drying column containing packed cotton wool, anhydrous sodium sulphate and silica. The extract was 
collected in the vial and concentrated by boiling it down with nitrogen gas to 1.0mL. The remaining extract was mixed with 1.0mL of 
the solvent and 1.0uL was injected into the flame ionization detector gas chromatograph for the TPH analysis. 
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Developing the Biokinetic Models 
The possible biochemical reaction to be observed as a result of micro-organism in a crude oil contaminated site is as shown below: 


Product oil + O2 + Nutrients ——-® CO>2 + H2O0 + New Biomass + Heat released 


From the equation shown above, the following mathematical models can be developed to illustrate the kinetics of production of 


biogenic gases produced in a reservoir for biodegradation of the soil. 


The Microbial Kinetics 
When a microorganism undergoes the biochemical reactions for dynamic studies, the general conservation for a steady state must 
be modified to give the following unsteady state mass balance. 


[nore _ [voworene _ rateof oil = a es (1) 


intothe outofthe : : : 
consumptionbymicrobes inthesystem 
system system 


For a batch system, the following assumptions hold true: 
1. The rate of oil flow into the system = 0 
2. Rate of oil flow out of the system = 0 


Hence equation (1) can be reduced to: 
rateofoil _ gens 
consumptionbymicrobes} inthesystem 


(2) 


ds 1 x aXe (3) 


ds 1 y Xt (4) 


ax _ ay, (5) 


Hence equation (4) becomes; 


aS _ __uXe 
dt Yx (6) 
5 


Monod's equation states that; 
m5: 
p= HS (7) 


Equation (6) can now be rewritten as; 
ds _ HmStXt (8) 


dt Y¥x(Km+S) 
s 


Where, x is the specific growth rate of microbes in the system per week, S, is the quantity of crude oil available at time t, X, is the 
microbial count at time t, Yy,s is the cell yield constant, < is the rate of depletion of the oil, < is rate at which the microbes 


proliferates, 41, is the maximum specific growth rate of the microbes in the system and Kj, is the half saturation constant 


Determining the values of 1, Kmand Yyx;s 
To determine 41, and Km, a graph of 1/u against 1/S was plotted. This implies that equation (7) had to be converted to an inverted 


form to achieve a straight line graph, i.e. 
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Lu Lm St Um 


Where, 


amis the slope of the graph and = is the intercept 


First, the values of were calculated using equation (10) below: 
n= n(8) x3 0) 
Where, X_ and Xare the initial microbial count and microbial count at time t respectively. 
Y,/s was determined using equation (11) below: 
X= Xs 


Ye/s = (11) 


So-St 


Make X; the subject of formula of equation (11), we have 
Xt = Ye/s(So — St) + Xo (12) 


Substituting equation (12) into equation (8), we have 


ds _ bmx /sSt—MmSo¥x/sSte~HmXoSt (13) 
dt YajsKm+¥a/sSt 


Evaluation of Equation (13) 
Equation (13) was solved using analytical methods to yield equation (14) below: 


a,log.(S — az) — az3log,S —a,—t =0 (14) 


Where, a1, a2, a3 and aq are arbitrary constants with values 30.9397599, 1371.64227752, 49.058803, -191.2084567 respectively for 
SRD and 19.6768957, 4567.1603183, 43.7196450 and -261.3550175 respectively for LRD. The value constants were obtained by using 
the collocation method and Newton-Raphson’s method. S is the Total Petroleum Hydrocarbon content. 


3. RESULTS AND DISCUSSION 


Data Analysis and Results 

This research work deals with the experimental results obtained from tests carried out, which was used in the computation to 
determine the functional parameters and coefficients. These results will lead to computations on the determination of maximum 
specific rate, specific rate and equilibrium constant. Results and data are represented on tables and figures. The figures will show 
similarity in the determination of maximum specific rate, specific rate and equilibrium constant for sandy soil and loamy soil each 
inoculated with 100g room-dried Azadirachta indica stem. Soil samples inoculated with 100g room-dried Azadirachta indica because 
they showed more promise than the rest set-up. The developed kinetics can be used to predict the entire biodegradation process of 
bonny light crude oil in Nigeria. 


Microbial Growth Rate (1) for SRD 

was calculated using equation (10), we have. 

For day 0 to day 7 where Xo is 3.94x 10cfu/g,Xt is 5.94x 10*cfu/g and t is 7, y is 0.374cfu/g/week 
For day 7 to day 14 where Xo is 5.42x104cfu/g,Xt is 1.03x 10’cfu/g and t is 7, p is 0.371cfu/g/week 
For day 14 to day 21 where Xo is 1.03 10’cfu/g,Xt is 1.99x 10%cfu/g and t is 7, p is 0.251cfu/g/week 
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For day 21 to day 28 where Xois 


1.99x 10°cfu/g,Xt is 5.99x 10'°cfu/g and t is 7, p is 0.12cfu/g/week 


For day 28 to day 35 where Xo is 5.99x 10'°cfu/g,Xt is 3.99x 10''cfu/g and t is 7, p is 0.05cfu/g/week. 


120 


100 


1/p 
fon) 
[o) 


1/uvs 1/S 


0.00E+00 2.00E-04 4.00E-04 6.00E-04 8.00E-04 1.00E-03 1.20E-03 


1/S 


Fig 1: Line Weaver-Burk Plot for SRD 


From Microsoft Excel, the equation of the plot in Figure 1 was found to be: 


7 = 89855.77~ + 24.038 


Comparing equation (15) to equation (9), we havey,, calculated to be 0.0416cfu/g/week and K,,to be 37380mg/kg. Equation 
(11) was used to obtain the value of Y,/;as2.596 x 10°cfu/g for X: = 1.987 x 10'4, Xo = 8.113 x 1013, So= 35818.69 and S; = 1003.18. 


(15) 


Having obtained the values of constantsy,,, Kmand Yx/s equation 12 can be represented as; 


ds _ 1.079x10°sS?—3.868x1017S;—3.375x 107? S; (16) 


dt 9.703X10124+2.596x109S; 


Microbial Growth Rate (1) for LRD 


1/p 


40 1/pu VS 1/S 


30 
20 
10 


0 
0.00E+00 5.00E-05 1.00E-04 1.50E-04 2.00E-04 2.50E-04 3.00E-04 
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Fig 2: Line Weaver-Burk Plot for LRD 


uy for LRD was calculated using equation (10), we have 
For day 0 to day 7 where Xo is 2.02 10'°cfu/g,Xt is 3.07x10'°cfu/g and t is 7, p is 0.0475cfu/g/week 
For day 7 to day 14 where Xo is 3.07x 10'°cfu/g,Xt is 4.22 x 10cfu/g and t is 7, u is 0.0455cfu/g/week 
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For day 14 to day 21 where Xo is 4.22x10'°cfu/g,Xt is 5.71 10"°cfu/g and t is 7, p is 0.043 1cfu/g/week 
For day 21 to day 28 where Xois 5.71x10'°cfu/g,Xt is 7.471x 10"°cfu/g and t is 7, p is 0.0384cfu/g/week 
For day 28 to day 35 where Xo is 7.471 x10'°cfu/g,Xt is 9.006x 10'cfu/g and t is 7, p is 0.0267cfu/g/week 


From Microsoft Excel, the equation of the plot in Figure 2 was found to be: 


: = 67291 ~ + 19.268 (17) 


Comparing equation 17 to equation 9,u,,, was calculated to be 0.0552/week and K,, to be 3713.8mg/kg. Equation (11) was used 
to obtain the value of Y,y/,as4.67 x 105cfu/g for X: = 9.006 x 10’, Xo = 2.202 x 101, So = 48508.69 and S; = 3707.78. Having 


obtained the values of constantsu,,, Kmand Ys equation 12 can be represented as; 


ds _ 2.58x104S?—9.23x10°S;—-1.214x10°S; (18) 
dt 1.73X109+4.67X105S; 


Comparison between the experimental TPH and Theoretical TPH 
MATLAB was able to estimate the trend of TPH depletion with time. The calculated values of TPH for 35 days are presented in Table 


2 below: 


Table 2: Comparison between Theoretical Valuesand Experimental Values of SRD and LRD 


: TPHih SRD TPHexp TPHexp LRD 
Time (days) TPHihLRD (mg/kg) 
(mg/kg) SRD(mg/kg) (mg/kg) 
35818.69 48508.92 35818.69 48508.92 
23477.06 35707.70 26865.7 37851.09 
14 15012.26 24890.51 14658.37 24890.09 
21 9085.47 16999.28 9085.48 16998.99 
28 3615.11 9875.66 3615.07 9897.95 
35 1003.33 3708.81 1003.18 3707.78 
TPH vs T 
60000 
50000 
B 40000 
~S 
0 
€ 30000 
ae 
a 
F 20000 
10000 
0 
0 5 10 15 20 25 30 35 40 
time (days) 
——TPHth SRD ———TPHth LRD (ppm) —-—= TPHexp SRD(ppm) ———= TPHexp LRD (ppm) 


Fig.3: Comparison between Theoretical Values and Experimental Values of SRD and LRD 


Using Analytical methods, equations (16) and (18) were solved to develop models that predict values of TPH of SRD and LRD 
respectively. Comparisons were made between theoretical values and experimental values using the percentage error equation: 
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TPHexp-TPHth 
TPHexp 


P.E= x 100(18) 


Where, P.E is the percentage error, TPHexpis the experimental value of TPH and TPHinis the theoretical value of TPH estimated by the 
new model 


The Table 3 below compares estimated values of TPH for SRD to their experimental equivalent. 


Table 3: Comparison of Experimental Values of TPH for SRD to the Calculated Values 


Time (days) TPHexp (mg/kg) TPHin (mg/kg) E (%) 
0 35818.69 35818.69 0 

7 26865.7 23477.06 12.6 
14 14658.37 15012.26 2.4 

21 9085.48 9085.47 0.004 
28 3615.07 3615.11 0.016 
35 1003.18 1003.33 0.014 


The Table 4 below compares estimated values of TPH for LRD to their experimental equivalent. 


Table 4: Comparison of Experimental Values of TPH for LRD to the Calculated Values 


Time (days) TPHexp (mg/kg) TPH cai (mg/kg) E; (%) 

0 48508.92 48508.92 0 

7 37851.62 35707.70 2.7 

14 24890.09 24890.51 0.00156 
21 16998.99 16999.28 0.00417 
28 9897.95 9875.66 0.00293 
35 3707.78 3708.81 0.000809 


For SRD, TPH at 21 days was calculated to be 9085.47mg/kg which is similar to that gotten from experiment observed to be 
9085.48mg/kg. The highest error was observed in day 7where calculated TPH was 23477.06mg/kg while the experimental value was 
26865.7mg/kg. 

For LRD, TPH at 35 days was calculated to be 3707.81mg/kg which is similar to that gotten from experiment observed to be 
3707.78mg/kg. The highest error was observed in day 7 where calculated TPH was 35707.07mg/kg while the experimental value was 
37851.62mg/kg. 


Predicting TPH during Bioremediation using the Developed Model 

Equation (14) in research work can be used to predict the TPH of the soil at any given time during bioremediation. The model was 
able to predict the TPH at Day 5, 10, 15, 30, 50, 55 and 60 picked arbitrarily. Table 5 and Figure 4 were used to illustrate the trend for 
both SRD and LRD. 


Table 5: TPH Prediction for Arbitrary Days 


Time (t) TPH-SRD (mg/kg) TPH-LRD (mg/kg) 
5 27538.03 38531.61 
10 19461.27 23644.18 
15 14037.14 18023.21 
30 1040.35 3845.07 
50 202.19 620.41 
55 187.20 589.11 
60 171.79 522.65 
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TPH vs Time 


TPH (mg/kg) 
N 
° 
fo} 
3S 
ro) 


-5000 0 10 20 30 40 50 60 70 
Time (days) 


—@—TPH-SRD —@=TPH-LRD 


Fig. 4: A Graph of TPH Predictions for Arbitrary Days 


The graph in Figure 4 established the accuracy of the developed predictive model. An attempt was made at interpolation using 
the predictive model, to find out TPH measures for arbitrary days within the period of experiment, thereby evaluating their 
accordance to the trend observed in the experiment. Sequel to day 7 of SRD whose TPH measure is 26865.7mg/kg from experiment, 
day 5 predicted by the model was seen to be 27538.03mg/kg indicating reasonable precedence. Also, between the theoretical TPH 
at day 5 and experimental TPH at day 7 for LRD, reasonable precedence was exhibited as the former has a TPH measure of 
38531.61mg/kg and the later has a TPH measure of 37851.09mg/kg. Attempt was also made at extrapolation using the predictive 
model, to check the consistency between TPH for arbitrary days outside the range of days used in the experiment and the TPH 
measure for the last day of the experiment. The model predicted that at day 50, the TPH measure for SRD would decline to 
202.16mg/kg which is a massive decline from the experimental 1003.18mg/kg at day 35. Also for LRD, the TPH measure predicted by 
the model would be 620.41mg/kg at day 50 which is a massive decline from the experimental 3707.78mg/kg at day 35. Hence, the 
model developed is reliable in predicting the trend of crude oil bioremediation using crushed room-dried Azadirachta indica stem. 


4. CONCLUSION 
A new model was developed to simulate the bioremediation of crude oil contaminated soils using Azadirachta indica stem. The 
experiment part of the study established the fact that 100g of blended room-dried Azadirachta indica stem caused the highest 
remediation in either sandy soil or loamy soil contaminated with 100m! of crude oil. Based on this observation, a biokinetic model 
was developed for both contaminated sandy soil and loamy soil treated with 100g blended room-dried Azadirachta indica stem. 
These parameters, together with Monod’s equation and material balance equation were then used to develop a model that 
simulates and predicts the process of crude oil bioremediation using Azadirachta indica stem. The TPH values obtained from the 
predictive model developed were similar to the TPH values gotten from the experiment. For SRD, TPH at 21 days was calculated to 
be 9085.47mg/kg which is similar to that gotten from experiment observed to be 9085.48mg/kg. The highest error was observed in 
day 7where calculated TPH was 23477.06mg/kg while the experimental value was 26865.7mg/kg. For LRD, TPH at 35 days was 
calculated to be 3707.81mg/kg which is similar to that gotten from experiment observed to be 3707.78mg/kg. The highest error was 
observed in day 7 where calculated TPH was 35707.07mg/kg while the experimental value was 37851.62mg/kg. 
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